Trial application of the envelope method to the potential ambiguity problem* 
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The optical potential ambiguity is a long-standing problem in the analysis of elastic scattering data. 

For a specific colliding system, ambiguous potential families can lead to different behaviors in the nearside 
and farside scattering components. 

By contrast, the envelope method can decompose the experimental data into two components with negative 
and positive deflection angles, respectively. Hence, a question arises as to whether the comparison between 
the calculated nearside (or farside) component and the derived positive-deflection-angle (or negative-deflection- 
angle) component can help analyze the potential ambiguity problem. In this study, we conducted a trial ap- 
plication of the envelope method to the potential ambiguity problem. The envelope method was improved by 
including uncertainties in the experimental data. 

The colliding systems of '©O+?8Si at 215.2 MeV and '*C+'°C at 1016 MeV were considered in the analyses. 
For each colliding system, the angular distribution experimental data were described nearly equally well by 
two potential sets, one of which is “surface transparent” and the other is refractive. The calculated angular 
distributions were decomposed into nearside and farside scattering components. Using the improved envelope 
method, the experimental data were decomposed into the positive-deflection-angle and negative-deflection-angle 
components, which were then compared with the calculated nearside and farside components. The capability of 


the envelope method to analyze the potential ambiguities was also discussed. 
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I. INTRODUCTION 


Optical potential represents the interaction between collid- 
ing nuclei pairs and plays an important role in the analyses 
of nuclear reactions [1]. The parameters of empirical optical 
potentials are typically obtained by fitting the experimental 
data of elastic scattering angular distributions (e.g., Refs. [2— 
5]). However, the ambiguity problem usually exists in the 
obtained potential parameters, which implies that the exper- 
imental data can be fitted equally well to different sets of 
potential parameters [1]. Potential ambiguities can be clas- 
sified into several different types, including continuous, dis- 
crete, refractive or diffractive, and shallow- or deep-W ones, 
which have been investigated in many studies (e.g., Refs. [6- 
10]). Clearly, the ambiguity introduces problems in analyzing 
the reaction mechanisms [11-15]. Therefore, proper treat- 
ment is required to resolve the ambiguity problem and ob- 
tain more physical optical potentials. One approach is to use 
global energy-dependent optical potentials (such as those in 
Refs. [16—-26]) to constrain the potential parameters. In ad- 
dition, for the elastic scattering of two light heavy-ion sys- 
tems, it was found that the refractive farside scattering data 
can help eliminate the discrete and shallow- or deep-W am- 
biguities [1]. 

In Ref. [10], Hussein and McVoy highlighted that at ener- 
gies well above the Coulomb barrier, the nearside and farside 


* This work was supported by the National Natural Science Foundation of 
China (Grant Nos. 12005047 and U1832105). 
t Corresponding authorsongyushou80@ 163.com 


scattering components exhibit different behaviors for the re- 
fractive and diffractive optical potentials. For the refractive 
potential, the nearside scattering component shows a larger 
slope than the farside component (when using In(do/d@) as 
the vertical coordinate). Consequently, the two components 
can construct the localized Fraunhofer crossover point, and 
a deep minimum appears near the corresponding scattering 
angle in the total differential cross sections. However, for 
the diffractive potential, the nearside and farside components 
are nearly parallel. In this case, rather than the localized 
Fraunhofer crossover point, a much wider Fraunhofer oscil- 
lation structure exists. The authors of Ref. [10] have taken 
the elastic scattering of 160+?8Si at 215.2 MeV as an ex- 
ample to show the differences between the nearside/farside 
decompositions calculated with the refractive potential and 
the diffractive one. For this colliding system, both the re- 
fractive deep-V potential called “A-type” potential and the 
diffractive shallow-V potential called “E18” potential (which 
is also known as “surface transparent”) can effectively repro- 
duce the existing experimental data of the angular distribu- 
tion [9]. However, the nearside and farside scattering com- 
ponents calculated using the two types of optical potentials 
exhibit different behaviors. For the nearside components, the 
difference in slopes is significantly small at angles smaller 
than those at which the nearside and farside components are 
close to each other. However, at larger angles, the nearside 
component calculated with the “E18” potential shows a rela- 
tively smaller slope than the one calculated with the “A-type” 
potential. For the farside scattering, different from the near- 
side case, the calculated cross sections with the “A-type” po- 
tential shows a smaller slope than that calculated with the 
“E18” potential in the whole angular range. These results 


indicate that a comparison between the nearside and/or far- 
side components can be used to demonstrate the differences 
between the ambiguous potential families. 

In Ref. [27], da Silveira and Leclercq-Willain proposed an 
envelope method to decompose the experimental data of the 
elastic scattering angular distribution into two scattering com- 
ponents with negative and positive deflection angles, respec- 
tively. In this method, two envelopes were drawn for the ex- 
perimental data with the well-defined maxima and minima in 
the data. Moreover, “experimental” data points correspond- 
ing to the negative- and positive-deflection-angle components 
can be derived. In Ref. [27], this method was applied to the 
elastic scattering data of a+4°Ca at 104 MeV and a+°8Ni at 
140 MeV without considering the uncertainties introduced by 
the experimental data. In addition, this method has been ex- 
tended to investigate the refractive phase relationship between 
elastic and inelastic scattering and the decomposition of the 
farside scattering angular distribution [28, 29]. In Ref. [30], 
the experimental dataset of a+*°Ca at 104 MeV was further 
analyzed using the envelope method. The obtained differen- 
tial cross-section data of the positive- and negative-deflection- 
angle components agree well with the nearside and farside 
scattering components, respectively. Combined with the be- 
haviors of nearside and farside scattering mentioned above, 
our observations inspired us to consider the envelope method 
as a potential solution to address the ambiguity problem. This 
could be achieved by comparing the “experimental” positive- 
and negative-deflection-angle components with the theoreti- 
cal nearside and farside components, respectively. 

In this study, a trial application of the envelope method to 
the potential ambiguity problem was performed. The exper- 
imental data of the elastic scattering angular distributions of 
160+?8Si at 215.2 MeV and 1?C+!?C at 1016 MeV were se- 
lected, and the optical potential ambiguity between the “sur- 
face transparent” and refractive potentials were studied. Op- 
tical model calculations and corresponding nearside/farside 
decompositions were performed for each dataset. Using the 
envelope method, the experimental data points within the an- 
gular range covered by the well-defined maxima and min- 
ima were decomposed into positive- and negative-deflection- 
angle cross-section components. To estimate the uncertainty 
of the two components, the errors of the experimental data 
points used in the construction of the envelopes were con- 
sidered. The capability of the envelope method to ana- 
lyze the potential ambiguity problem was discussed by com- 
paring the obtained positive/negative-deflection-angle cross- 
sectional components with the calculated results of the near- 
side/farside decompositions. 


Il. METHOD 


The optical model, the nearside/farside decomposition and 
the envelope methods were utilized to analyze the experimen- 
tal data. In this section, the latter two methods are briefly 
outlined; more details can be found in Ref. [27, 30-33]. 

In the envelope method, the starting point is a partial-wave 
expansion of the scattering amplitude [27]. Using the asymp- 


totic expression of the Legendre polynomials, the scattering 
amplitude can be decomposed into the components corre- 
sponding to the positive and negative deflection angles, which 
are denoted by f,(0) and f_ (6), respectively. 


F) = f4(0) + f- (0). (1) 


In the classical limit, f,;(0) and f_(@) can be expressed as 
follows [30]: 
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(2) 
where ô; is the nuclear plus the Coulomb phase shift. The 
cross sections corresponding to the positive and negative de- 
flection angles are denoted as o4 (0) and o— (0), respectively, 
which are related to the scattering amplitudes as follows: 


o+ (8) = |f+(9)/?. (3) 


It should be noted that the signs (“+” and “—”) used as sub- 
script are similar to those used in Ref. [30] but in contrast to 
those used in Ref. [27]. The elastic scattering angular dis- 
tribution o()=| f(0)|? oscillates between the upper envelope 
E, and lower envelope EF; defined by 


r 


E, (0) = (a 2(@) 4 HO l (4) 
and 


2 
E,(6) = (030) So °(6)) (5) 
According to Eq. (4) and (5), the positive- and negative- 
deflection-angle scattering cross sections can be obtained us- 
ing two envelopes of the experimental elastic scattering an- 
gular distribution. The upper and lower envelopes are deter- 
mined by the maxima and minima of the angular distribution, 
respectively. The envelopes can be drawn if the maxima and 
minima are well-defined in the experimental angular distri- 
bution. In the present study, the envelopes were drawn by 
connecting these maxima and minima with broken lines. The 
broken lines were then smoothed using the TGraph class pro- 
vided by the ROOT software [34]. 

Observing Eq. (4) and (5), one can find that c4 (0) and 
a—(0) follow the same relationship. Additional analysis is re- 
quired to determine whether the solution obtained is o4 (0) or 
a—(6@). This problem can be resolved by observing the trends 
of the ratio of o4 (0) to the Rutherford scattering differential 
cross section op(@) because o4 (0)/opr(0) tends to unity as 
0 approaches zero and decreases almost exponentially as 0 
becomes larger [27, 30]. 

In previous studies based on the envelope method, uncer- 
tainties in the experimental data have not been considered 
quantitatively. For each maxima/minima data point, the mean 
value (u) was used to determine the envelope. In the present 
work, in the construction of the envelopes, the standard de- 
viation (dev) of each maxima/minima data point was also 
considered. In addition to u, the envelopes were drawn for 


u + dev and u — dev to estimate their uncertainties. Hence, 
there were three upper and three lower envelopes. The mean 
values of o4 (0) were obtained by combining the upper and 
lower u envelopes. The other combinations of upper and 
lower envelopes were used to limit 74(@) and o_(@). Us- 
ing this approach, the uncertainties can be propagated to the 
calculated o+(0), which makes the envelope method more 
quantitative. 

The nearside/farside decomposition of the scattering am- 
plitude was performed using the method proposed by 
Fuller [31]. The scattering amplitude can be expressed as 


F(9) = fr(9) + Fr). (6) 


where fn(0) is the nearside scattering amplitude: 


fin (8) a 20 + 1)[e?** E 2% Q (cos) 
=0 


(7) 
+ fon(0), 
and f(@) is the farside scattering amplitude: 
1S l 
FO) = X2 + Ve" — e] (cos0) 
2ik - i (8) 


+ fo,+ (8). 


In Eq.(7) and Eq. (8), fc,n and fc,¢ are the nearside and far- 
side Coulomb scattering amplitudes, respectively; k is the rel- 
ative wave number; a; is the Coulomb phase shift; and the 


travelling wave functions Q (cos@) are expressed as 


Q (cos0) = S [Pi(cos0) + {Or(cos0) 


where P;(cos 0) and Q;(cos 0) are the Legendre functions of 
the first and second types of degree l, respectively. The differ- 
ential cross sections corresponding to the nearside and farside 
scattering amplitudes are denoted as o» and oy, respectively: 


on/p (9) = | fns? (10) 


In Ref. [30], it was shown that when Q is replaced with 
asymptotic forms, f, and fp are simply amplitudes f, and 
f_, respectively. 

In the present work, the comparisons between “experimen- 
tal” o+ (0) and the theoretical o, (@) were performed to test 
the capability of the envelope method in treating the poten- 
tial ambiguity problems. The optical model calculations were 
performed using FRESCO [35]. The nearside/farside decom- 
position was conducted according to the method described in 
Ref. [32]. 


HI. RESULTS AND DISCUSSIONS 


Experimental data for the elastic scattering angular distri- 
butions of 1©O+28Si at 215.2 MeV (Ref. [36]) and !2C+!C at 


TABLE 1. The nuclear optical potential parameters used in the op- 
tical model and nearside/farside decomposition calculations. The 
parameters for '©O+?°Si at 215.2 MeV and '?C+!?C at 1016 MeV 
are extracted from Ref. [9] and Ref. [11], respectively. 
colliding potential Vo TR aR Wo TT ay 
system set (MeV) (fm) (fm) (MeV) (fm) (fm) 
1694285; Al 10.00 1.350 0.618 23.40 1.230 0.552 
A2 100.00 0.967 0.745 44.10 1.073 0.850 
Bl 49.90 0.934 0.742 150.40 0.262 1.201 
B2 129.4 0.681 0.913 47.90 0.918 0.622 
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Fig. 1. Elastic scattering angular distribution of '©O+?8Si at 215.2 
MeV. The circles represent the experimental data. The curves are 
the results of the optical model and nearside/farside decomposition 
calculations with different potentials. 


1016 MeV (Ref. [37]) were analyzed using the optical model, 
the nearside/farside decomposition, and the envelope method, 
respectively. 

Experimental data for 1°O+?°Si at 215.2 MeV were ob- 
tained by digitizing the data points shown in “Fig. 6” 
in Ref. [38]. The digitizing software used in this study 
was GSYS (version 2.4.9) [39]. The experimental data 
for 1*C+!7C at 1016 MeV were obtained from the NRV 
database [40]. 

The optical potential parameters used in the optical model 
and the nearside/farside decomposition calculations are listed 
in Table 1. For both the systems, the “standard” Woods-Saxon 
function was adopted for the nuclear optical potentials: 


Vo i Wo 
1+ exp (28) 1 + exp (=) (11) 


U(r) = 
where Vo and Wo are the depths of the real (R) and imaginary 
(D) parts, respectively. R; and a; (i = R, J) are the radius and 


diffuseness parameters, respectively. The reduced radius r; 
(i=R, I) is defined as 


r; = R,/(Af? + AH’), i=R,I (12) 
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Fig. 2. Elastic scattering of '©O+?8Si at 215.2 MeV: (a) envelopes of the experimental data; (b) greater/smaller components in comparison 


with the nearside/farside components. 


where Ar and Ap are the atomic masses of the target and pro- 
jectile, respectively. For the Coulomb potential, the potential 
between a point charge and a uniformly charged sphere was 
utilized. The optical potential parameters for '©O+?°Si and 
12C+!2C were obtained from Ref. [9] and Ref. [11], respec- 
tively. 


A. 1°Q+?5Si at 215.2 MeV 


First, we revisited the refractive or diffractive potential am- 
biguity problem of 1°O+78Si at 215.2 MeV, as discussed in 
Ref. [9, 10]. In fact, the potential sets Al and A2 are simply 
the diffractive shallow-V “E18” and refractive deep-V “A- 
type” ones, respectively [9]. The A1 potential is also a “sur- 
face transparent” potential, whose characteristic is a larger 
value of Rpr than that of Ry. Its real part is deeper than 
the imaginary part for r > Rp [10]. The results of the op- 
tical model and nearside/farside decomposition calculations 
are shown in Fig. 1. It can be seen that both the potential 
sets Al and A2 can effectively reproduce the experimental 
data, although the theoretical angular distributions show quite 
different patterns at angles beyond approximately 35° where 
no data exist. As reported in Ref. [10], within the angular 
range covered by the experimental data, although the near- 
side components corresponding to Al and A2 potentials are 
significantly close to each other, the farside components are 
evidently different. 

The upper and lower envelopes of the experimental data 
are presented in Fig. 2(a). For each envelope, the uncertain- 
ties are represented by the region delimited by dotted lines. 
The positive- and the negative-deflection-angle components, 
o+(0) and o- (0), are derived by using the envelopes and are 
shown in Fig. 2(b). The envelope uncertainties are propagated 
to o4 (0). The uncertainties in o4 (0) are also represented by 
the regions delimited by dotted lines. For comparison, the 
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Fig. 3. Elastic scattering of *°C+*?°C at 1016 MeV. The circles 


represent the experimental data. The curves are the results of the 
optical model and nearside/farside decomposition calculations with 
different potentials. 


nearside and farside angular distributions, o,,(@) and a (@), 
are also plotted. It can be observed that 7+ (0) matches on (0) 
well with small uncertainties. Moreover, a_ (8) is close to 
the magnitude of o;(@). As the calculated on (0) with poten- 
tials Al and A2 are significantly close to each other, it is not 
feasible to determine which is closer to the derived o4 (0). 
However, the calculated o/(@) with potentials Al and A2 
are evidently different. However, the relatively small values 
(less than 107°) of o— (0) probably make its apparent struc- 
ture meaningless [28]. In addition, the relative error of o— (6) 
is significantly larger than that of o4 (0), which is consistent 
with the results calculated in Ref. [30], in which the author 
found that slight adjustments of envelopes visibly affect the 
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Fig. 4. Elastic scattering of '*C+!?C at 1016 MeV: a) envelopes of the experimental data; (b) greater/smaller components in comparison 


with the nearside/farside components. 


smaller ones of o4 (0) and o_ (0) but have little effect on the 
larger one. These difficulties render it impractical to deter- 
mine whether o f(@) is closer to o_ (0). Hence, the envelope 
method is unsuitable for resolving the refractive or diffractive 
potential ambiguity problem of !©O+78Si at 215.2 MeV. 


B. 1C+'?C at 1016 MeV 


For the elastic scattering of 17C+!2C at 1016 MeV, the cal- 
culated total, nearside and farside angular distributions with 
the B1 and B2 potentials are plotted in Fig. 3 in comparison 
with the experimental data. For the B1 potential, although 
the depth of the real part is smaller than that of the imagi- 
nary part, it is not a diffractive potential because it produces 
a Fraunhofer crossover [10]. Similar to the Al potential, the 
B1 potential is a “surface transparent” one with Rp > Ry. 
The strength of the real part is greater than that of the imag- 
inary part for r > 2 fm, and there is a strongly absorbing 
core at small radii [11]. The B2 potential is a refractive deep- 
V potential, similar to the A2 potential. In these calculations, 
following the procedure described in Ref. [11], the relativistic 
effects were considered by using an effective center-of-mass 
energy of 496.84 MeV and an effective mass of 12.27 u for 
12C, Observing Fig. 3, both the B1 and B2 potentials effec- 
tively reproduce the experimental data. The farside compo- 
nents g ¢(@) are significantly close to each other in the angular 
range of the experimental data. However, the nearside com- 
ponents o,,(@) exhibit an evident difference in slope after ap- 
proximately 5° in the center-of-mass system, in which o,,(0) 
corresponding to the B1 potential becomes smaller than that 
corresponding to B2. 

The upper and lower envelopes are drawn for the experi- 
mental data in Fig. 4(a), and the corresponding derived c+. (0) 
are drawn in Fig. 4(b) in comparison with o,,/7(0). As in 
the case of !©O+?8Si, the envelope uncertainties are also con- 


sidered and propagated into o(@). It can be observed that 
o_(@) matches o/(@), and o+(@) is close to o,,(@) for both 
B1 and B2 potentials. In this case, because the farside scat- 
tering components o (8) corresponding to the two potential 
sets are close to each other, it is not feasible to determine 
which one is closer to the derived o_ (0). However, the ev- 
ident difference between the calculated o,,(0) with the two 
potential sets enables the selection of the physical potential. 
As for the obtained c+ (0), its magnitude is larger than 107° 
at almost all the scattering angles covered by the envelopes. 
In addition, the relative error increases as the scattering angle 
increases. Clearly, o+ (0) matches better with the B2 nearside 
component than with the B1 component. This reveals that the 
B2 potential is more physical, which is consistent with the 
results of the combined analyses of elastic scattering and the 
transfer reaction of !2C+!2C at 1016 MeV in Ref. [11]. 


C. Factors influencing the application of the envelope method 


The calculated results for 1*C+!2C at 1016 MeV prelimi- 
narily indicate the possibility of using the envelope method to 
address the potential ambiguity problem. Several factors in- 
fluencing the application of the envelope method to potential 
ambiguity problems were identified. 

The primary factor is the difference between the calculated 
On/ (9) corresponding to different potential families. If they 
show evident differences in o,,(@) and/or of(@), there is a 
possibility of selecting a physical potential. However, when 
the calculated o,,/¢(@) values are nearly identical for different 
potential families, applying this method is not feasible. 

The second factor is the magnitude of the smaller part of 
o+(0). When the calculated o,,/; (0) exhibits evident differ- 
ences from the smaller values, as in the cases of !°O+?8Si 
at 251.2 MeV and 1?C+!?C at 1016 MeV, the magnitudes 
should not be significantly small. The present calculations 


show that when the magnitudes are smaller than approx- 
imately 1073, a comparison between the smaller parts of 
o+(@) and o,/(@) can provide limited information. In 
this case, although the derived o+ (0) matches the calculated 
On/f(@) well on the order of magnitude, the details of the 
structure become meaningless. This phenomenon was ob- 
served in the positive-/negative-deflection-angle decomposi- 
tion of the inelastic scattering angular distributions using the 
envelope method in Ref. [28]. 


Furthermore, the data quality also influences the applica- 
tion of the envelope method. The plot of the envelopes re- 
quires well-defined maxima and minima; in other words, it 
requires good angular resolution. The errors of the differen- 
tial cross sections should be small as they are propagated to 
the derived o4 (0). 

Apart from the factors encountered and mentioned above, 
it should be noted that the Fraunhofer crossover can also in- 
fluence the application. In fact, at the Fraunhofer crossover 
angle, the values of o (0) and o— (0) should equal with each 
other in principle, which implies that the value of the lower 
envelope should be zero (cf. Eq. (5)). However, the lower 
envelope does not have zero points. Therefore, as indicated in 
Ref. [30], the data near the Fraunhofer crossover are unsuit- 
able for use in the envelope method. 


IV. SUMMARY 


In this study, a trial application of the envelope method to 
the optical potential ambiguity problem is presented. The col- 
liding systems of 1©O+?8Si at 215.2 MeV and 1?C+!?C at 
1016 MeV were used as examples in the application. Based 
on the envelope method, the experimental data of the elas- 
tic scattering angular distributions were decomposed into 
positive- and negative-deflection-angle components (denoted 
as o4 (0) and o_(6), respectively). Although the original ver- 
sion of the envelope method does not quantitatively consider 
uncertainties, the present work propagates the uncertainties 
of the experimental data points that determine the envelopes 
into the derived o+ (0). For each colliding system, the op- 
tical model and the nearside/farside decomposition calcula- 
tions are performed with two potential sets, of which one is 
“surface tranparent” and another is refractive. For each po- 
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tential set, the corresponding nearside and farside scattering 
components (denoted as o,,() and (6), respectively) were 
compared with the derived o4 (0) and o_ (6), respectively. 

For 1©0+?8Si at 215.2 MeV, the envelope method provides 
good estimations of o,,/¢(@) on the order of magnitude. The 
nearside components o,,(@) corresponding to two potential 
sets are nearly identical in the angular region of the envelopes. 
Therefore, although the derived o4 (0) matches them well, 
it is not feasible to select a physical potential by compar- 
ing 74(@) with o,,(@). By contrast, the farside components 
a s (0) corresponding to two potential sets, exhibit evident dif- 
ferences in the angular region of the envelopes. However, the 
derived o_(@) has small magnitudes with large uncertainties, 
making its structure meaningless. It is impractical to select 
the physical potential by comparing o_ (0) with o¢(@) quan- 
titatively. 

For 12C+!°C at 1016 MeV, the derived o4 (0) and o_ (0) 
match the calculated o,,(@) and o/(@) well, respectively. In 
this case, the farside components o f(@) corresponding to two 
potential sets are significantly close to each other in the an- 
gular region of the envelopes. Although the derived o_ (0) 
matches well with a; (6), it is not feasible to select the phys- 
ical potential by comparing o—(0) with o/(6). For on(0), 
the two potential sets yield significantly different results. The 
calculated o,,(9) corresponding to the “surface transparent” 
potential B1 shows an evident smaller values after 5° in com- 
parison with that corresponding to the refractive potential B2. 
The derived o+(6) clearly matches much better with on (0) 
corresponding to B2. Hence, the refractive potential B2 was 
selected as the physical potential, which was consistent with 
the results given in Ref. [11]. 

Several factors that influence the application of the enve- 
lope method to the potential ambiguity problem were identi- 
fied. We believe that one can try to use the envelope method 
to analyze the potential ambiguity problem when the follow- 
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